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We calculate pseudoscalar masses on gauge configurations containing the effects of 2+1 flavors 
of dynamical asqtad quarks and quenched electromagnetism. The lattice spacings vary from 0.12 
to 0.06 fm. The masses are fit with staggered chiral perturbation theory including NLO electro- 
magnetic terms. We attempt to extract the fit parameters for the electromagnetic contributions, 
while taking into account the finite volume effects, and extrapolate them to the physical limit. 
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1. Introduction 



Quark masses are fundamental parameters in the Standard Model and their determination from 
lattice field theory is important for phenomenology. With the advance of high-precision spec- 
troscopy the electromagnetic (EM) contributions to the quark masses have become important. In 
fact, they constitute the largest uncertainty in the calculation of m u /md [|I|] as indicated in Table |]. 
The EM error stems from uncertainties in the size of the EM contributions to the masses of the 
71 and especially the K mesons. The contributions have until recently been taken from a vari- 
ety of phenomenological estimates, and therefore have quite large and not well controlled errors. 
The MILC collaboration has been working on reducing these uncertainties and progress has been 
reported previously in Refs. [§, HJ, |]]. 

To lowest order (LO) in chiral perturbation theory {%PT) the electromagnetic splitting in the 
pion and kaon systems is identical (this observation is known as Dashen's theorem [||]). We aim to 
calculate on the lattice the corrections to Dashen's theorem by employing staggered %PT with EM 
effects for fits to spectrum data. The corrections to Dashen's theorem may be parametrized as [a]: 



(M 2 K± -M 2 K0 y ={\ + e){M 2 n± 



(i.i) 



where the superscript / denotes the EM contributions. 1 Our main motivation for computing e is to 
be able to remove the EM effects from the experimental K masses; the resulting pure-QCD masses 
may be compared with QCD simulations in order to extract m u /md. Staggered %P1 ', developed 
in Ref. gives for the EM splitting between Goldstone pseudoscalars with charged and neutral 
quarks, AM 2 5 , the following NLO expression: 
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where a is the sea-quark label, £ is the staggered taste, q^ is the charge of a meson made of a 
quark x and an antiquark y (for more details on the notations see Ref. To fit our spectrum data, 



m u [MeV] m d [MeV] m u /m d 



value 


1.9 


4.6 


0.42 


statistics 


0.0 


0.0 


0.00 


lattice systematics 


0.1 


0.2 


0.01 


perturbative 


0.1 


0.2 




EM 


0.1 


0.1 


0.04 



Table 1: Value of m M /m^ calculated by the MILC Collaboration [[]]] and estimates of the different sources 
of error (with m u and mj determined in the MS scheme at 2 GeV). The row labeled "EM" shows the size of 
the error coming from a phenomenological estimate of electromagnetic contributions to the quark masses. 



We do not include quark-mass renormalization effects due to the presence of the EM field, since they are much 
smaller than our precision. 
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~ a [tin J 


Volume 


P 


mi j m s 


# configs. 


0.12 


20 3 x 64 


6.76 


01/0 OS 

\J . \J J. / \J . \J-J 


2254 




28 3 x 64 


6.76 


0.01/0.05 


274 




20 3 x 64 


6.76 


0.007/0.05 


1261 




24 3 x 64 


6.76 


0.005/0.05 


2099 


0.09 


28 3 x 96 


7.09 


0.0062/0.031 


1930 




40 3 x 96 


7.08 


0.0031/0.031 


1015 


0.06 


48 3 x 144 


7.47 


0.0036/0.018 


670 



Table 2: Parameters of the asqtad ensembles used in this study. 



we use the above expression plus the NNLO analytic terms. For investigating the finite volume 
effects, we include standard terms dependent on m n L s from EM tadpoles, as well as an empirical 
EM finite volume correction in the form f v q\ y /L 2 s used previously in Ref. [0], where f v is a constant 
and L s is the spatial lattice size. However, as discussed below, our measured finite volume effects 
are rather small at present, and including or omitting the finite volume terms from the fits currently 
makes little difference in the final results. 

2. Lattice setup 

Our EM field on the lattice is "photon-quenched," i.e., the sea quarks are taken to be electro- 
magnetically neutral, and the C/(1)bm links are generated independently of the SU{3)qcd links. 
The fact that the photons are not present in the sea does not affect the NLO fits to the EM meson 
splittings since, as shown in Ref. [||], the unknown LEC controlling sea-quark mass dependence 
cancels in EM splittings between mesons with the same valence masses. There are physical NLO 
effects missing due to the absence of sea-quark charges, but these may be inserted after the fact 
since they appear only in chiral logarithm terms with known (LO) coefficients. Uncontrolled EM 
quenching effects do come in at NNLO, but should be negligible at the current level of precision. 
Because diagrams that have intermediate gluon or photon states ("quark-disconnected" diagrams) 
are difficult to compute on the lattice, and because of the photon quenching, we have to redefine 
the neutral pion on the lattice '71°' to have the mass: 

= ( m la + m L")/ 2 ' ( 2 - 

where uu and dd are states with charges (2/3e, —2/3e) and (l/3e, —\/3e) and disconnected dia- 
grams are dropped by fiat. The effect of the disconnected diagrams is expected to be small, since 
EM contributions to the neutral pion must in any case vanish in the chiral limit. 

We calculate the meson spectrum on the set of asqtad ensembles with 2+1 flavors listed in 
Table | The values of the electron charge that we use are e = ±0.909, ±0.606 and,±0.303,0, 
corresponding to 9a p h ys , 4a p h ys , cc p hys an( l 0- 

3. Results 

We do not see large finite volume effects in the EM splittings calculated on the a « 0. 12 fm 
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ensembles with L s = 20 and 28. The volume effects do increase with the quark charges, but remain 
small. In Fig. [j] we compare our results with the finite volume effects for a charged kaon (with 3 
times the physical charges) expected from the BMW collaboration result [0] (magenta error bar). 
We find that ours are about three times smaller. A similar comparison is shown in Fig. |3| for 
physically charged kaons. Nevertheless, the error on the finite volume effect in our data (computed 
simply by propagating the errors in the difference between the results on the two volumes) is 
slightly larger than the effect itself, so any tension with the BMW result is less than a 2a effect. 
Indeed, we see no reason for any difference with BMW other than statistics, since as far as we can 
tell, we are treating the £/(1)em fields in the same way they do, and we do not expect the EM finite 
volume effects to be strongly dependent on quark mass or quark action. The finite volume effects 
we see are also smaller than those expected from the RBC Collaboration results [Q], but in that case 
some or all of the difference may be due to the different treatment of zero modes in our £/(1)em 
fields. We are currently increasing the statistics on the a w 0.12 fm, L s = 28 ensemble in hopes of 
clarifying this issue. 

The taste splitting between the Goldstone pion and the second local pion grows with the quark 
charge for the neutral pions (see Fig. || for the two coarse ensembles with the same parameters and 
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Figure 1: EM pseudoscalar splittings calculated on two coarse ensembles with mi/m s = 0.01/0.05 and 
spatial sizes L s — 20 and L s = 28 vs. the sum of the valence quark masses. The results for the different quark 
charge combinations are distinguished by color. The size of the finite volume effects found by the BMW 
group [0] is shown in magenta in the upper right and is scaled to charge 3e, to correspond to our magenta 
points plotted slightly below that. The finite volume effects we see (difference between crosses and squares) 
are smaller than those found by BMW. 
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Figure 2: Taste splitting between the Goldstone pion and the second local pion on the coarse ensembles 
with in/ jm s = 0.01/0.05 and L s = 20 and L s = 28. The x-axis values are arbitrary labels for the mesons with 
various values of the quark charges. Case 1 is a pion made from uncharged quarks. The first group of neutral 
pions (cases 2-4) are essentially uu combinations, and the second neutral group (cases 8-10) are of the dd 
type (the (2e/3, — 2e/3) combination appears in both neutral groups). These groups are labeled by one of 
the quark charges (the other has an opposite sign). The middle group (cases 5-7), has quark charges with 
the physical ratio of (2: 1); they are labeled by the total meson charge. 



different L s in Table |2|). This shows that, for charges larger than the physical charges, quark taste 
changing by high-momentum photons becomes a significant correction to that by high-momentum 
gluons. For charged pions, where annihilation into photons is forbidden, the effect is smaller (and 
negative), but still evident. The fit function Eq. ( |L2| ) is based on the neglect of taste-violations 
caused by photons. For that reason, the analysis presented here focuses only on fits with physical 
quark charges. However, to the extent that photon-induced taste violations remain small compared 
with gluon-induced violations, which is arguably the case for all values of quark charges shown 
in Fig. ^, one should be able to expand the fit function in powers of (Xem- This would allow 
such effects to be included by adding a| M analytic terms to the fit function. We have made some 
preliminary fits along these lines to our full data set with quark charges up to three times the 
physical values, and they seem to work reasonably well. 

Before leaving the subject of taste violations, we note that Fig. ||] shows that taste splitting is 
not very much affected by finite volume effects (compare the data in blue and red at L s = 20 and 
28 for the above coarse ensembles). 

We now fit, to Eq. ( |L~2| ) plus NNLO analytic terms, the EM splitting of the meson spectrum 
using the data with the physical electron charge e (i.e., mesons made from combinations of u 
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and/or d quarks, with q u = 2/3e and qj = — l/3e). Partially quenched points are included; in 
typical fits there are between 50 to 120 data points, with 20 to 30 fit parameters (depending on how 
many of the NNLO terms are included, and whether small variations with a 2 of the LO and NLO 
low energy constants are allowed). A typical fit is shown in Fig. ||, along with the unitary pions 
and kaons on each ensemble. The fit is uncorrected: the covariance matrix is nearly singular, 
and the statistics are insufficient to determine it with enough precision to yield good correlated 
fits. Based on the fits, our preliminary determination of the correction to Dashen's theorem in the 
continuum limit and at physical quark mass is e = 0.65(7)(14), where the errors are statistical 
and systematic, respectively. The systematic error comes largely from variations in e when the 
assumptions going into the chiral/continuum fit are changed. Systematic effects are still under 
investigation. In particular, based on the currently small finite volume effects seen in our data, 
we do not include a finite volume error here. However, we suspect that better statistics on our 
larger lattices will yield a statistically significant finite volume error. Our result is compatible with 
previous ones: 0.60(14) (statistics only) @], 0.628(59) (statistics only) [g], and 0.70(4)(8) [@]. 
With our value for e from above, our preliminary estimate for the EM uncertainty in m u /md is 
reduced to 0.022 [11 1, which is about half of our previous error shown in Table |]. 
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